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ABSTRACT 

The highly organised distribution of satellite galaxies surrounding the Milky Way is a serious 
challenge to the concordance cosmological model. Perhaps the only remaining solution, in 
this framework, is that the dwarf satellite galaxies fall into the Milky Way's potential along 
one or two filaments, which may or may not plausibly reproduce the observed distribution. 
Here we test this scenario by making use of the proper motions of the Fornax, Sculptor, Ursa 
Minor and Carina dwarf spheroidals, and trace their orbits back through several variations 
of the Milky Way's potential and account for dynamical friction. The key parameters are 
the proper motions and total masses of the dwarf galaxies. Using a simple model we find 
no tenable set of parameters that can allow Fornax to be consistent with filamentary infall, 
mainly because the la error on its proper motion is relatively small. The other three must 
walk a tightrope between requiring a small pericentre (less than 20 kpc) to lose enough orbital 
energy to dynamical friction and avoiding being tidally disrupted. We then employed a more 
realistic model with host halo mass accretion, and found that the four dwarf galaxies must 
have fallen in at least 5 Gyrs ago. This time interval is longer than organised distribution is 
expected to last before being erased by the randomisation of the satellite orbits. 

Key words: galaxies: dwarf galaxies: evolution gravitation Local Group dark matter 
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1 INTRODUCTION 

The observed number, spatial distribution and orbital energy 
of satellite galaxies surrounding the Milky Way is quite surpris- 

■ ing in terms of th e concordance cosmological model (see e.g. 

' lSDergeletal]|2007h . where galaxy hosting cold dark matter ha- 
los are formed by the hierarchical merging of many subhalos. Be- 
fore any baryonic physics has been considered, we would expect 
a moderately anisotropic distribution of cold dark matter subha- 
los surrounding the main halo, like shown in 'Li & Helmi' ('20081), 
that number in the hundreds (see in particular Klvpin et al. 1999; 
lMooreetal.1 19991) . This prediction is clearly at odds with the small 
~ 0(10) number of observed satellites - which shou ld no t include 
the unbound Magellanic C lou ds: see e.g.[K allivavali l et al. H2006ah . 
iKallivavalil etaU ( l2006bl) and lPiatek et all i2007bh . 

The perfectly mundane solution to this problem is that 
the majority of the satellites are not efficient at forming stars 
and so stay dark. This is more or less achievable if one tunes 
the available free parameters surrounding the feedback of su- 
pern ovae and other poorly understood astrophysical processes 
(lEfstathiou,200Q : .Madau et aL,2001, : .Springel & Hernquist.2003. : 



E-mail: angus.gz@gmail.com 



IVeilleux et ai]|2005l : loovernato et al.ll2007l amongst many others) 
contributing to star and galaxy formation. 

This is not the end of the story, however, because the spa- 
tial distribution of these remai ning satellit es is equally puzz ling. 
As highlighted in iKroupaetal] 12 005): Me tz & Kroupal ( l2007h and 
previously noted by iLvnden-Belll ( 119831) , their current positions 
around the Milky Way can be likened to an extended disk with a 
root mean squared thickness of merely 10 — 30 kpc. This is again 
contrary to the naive expectation in the concordance cosmological 
model because a priori one would expect the luminous satellites to 
be randomly distributed. 

This argument h as been countered bv lLi&Helmil ( l2008h (see 
also lD'Onghiall2008h who claimed the thin nature of the disk of 
satellites (DoS) might be explained by the infall of a small group of 
galaxies that wou ld retain the correlated orbi ts. T his suggestion was 
dissected by both Klimento wski et al.l ^20 id) and Met z et"ai] j2009l) 
who both gave arguments as to why the model would likely fail. In 
particular. iMetz et"^ ^20091) pointed out that no nearby groups are 
obs erved to be anywher e near as flattened as the disk of satellites, 
and lKlimentowski et al. showed, with N-body simulations, 

that the coherence timescale for these groups is short lived after 
infall. 

It is worth keeping in mind, however, that the Milky Way 
is not isolated and that the presence of M3 1 in our vicinity may 
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have some beari ng on the formation of the Milky Way system a s 
a whole (see e.g. lSawa & Fuiimotolliool : IPasetto & Chiosill2007h . 
Also there is the possibility of an intricate set of circumstances 
whereby several dwarfs in a dynamical group are sacr ificed to al- 
low th e survival of single member through scattering jSales et al] 
l2007al lbl). Howe ver, perhaps besides these, the only remain i ng so- 
lution is that oflKangetalJ feOOSh : iLibeskind et"5] j2005l. I200I 
l2009l) : IZentner et alj feooj ) who have suggested that the satellites 
fall in at high redshift along filaments which creates a preferred or- 
bital direction and might satisfy the small width of the DoS. How- 
ever, there is a further piece of the jigsaw that has been neglected, 
or overlooked until now, and that is the proper motions of the dwarf 
satellites. When these dwarfs fall along the filament, they accelerate 
towards the Milky Way and, to match the observed proper motions, 
they must lose orbital energy along the way, presumably to dynam- 
ical friction. 

This is important because the proper motions of four of the 
classical dwarf spheroidal galaxies (dSphs) - Forn ax, Sculptor, 
Ursa Minor and C a rina - have been measured by IPiatek et"ail 
\2miA I2OO6I l2005h ; iMetz et al.l ( l2008b respectively. In order for 
the filamentary infall hypothesis to be valid it must be possible for 
satellites to fall in to some realistic Milky Way potential and at 
some point reach the current positions of each of the dSphs and be 
consistent with their measured proper motions and radial velocities. 
The purpose of this paper is to test this argument. 



2 INITIAL CONDITIONS AND STRATEGY 

The strategy we employ to test the filamentary infall scenario 
is quite simple. 



2.1 Dynamical Friction 

We modified an orbit integrator, which usually simply updates 
the position and velocity of a test particle according to some grav- 
itational potential, to include the effects of dynamical friction. Dy- 
namical friction (DP) is the retarding force experienced by a mas- 
sive object as it orbits through a sea of particles or stars. As it 
moves, the object scatters oncoming particles creating an overdense 
wake to its rear. This overdensity acts gravitationally to decel erate 
the massive object. According to Isinnev & Tremaind bOOSh . the 
deceleration due to DF can be parametrised as follows 
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Here, Md is the total mass of the dSph and V is the velocity of 
the dSph with respect to the Milky Way; G is Newton's constant, 
Vc (r) is the local circular speed defined by the Milky Way potential, 
Pdm{r) is the local dark matter density (we assume that the dSphs 
orbit outside the disk) and ln(l -|- A^) is the Coulomb Logarithm. 
Ordinarily, ln(l + A^) can be reduced to 2 In A, but in our case A 
can be quite small, so it is better left unaltered. 

In addition to this, A is often parametrised as ''gg^^ where 
bmax is linked to the maximum impact parameter of particles 
that meaningfully contribute to the deceleration. In this work we 
take bmax ~ 20 kpc, which is 20 — 40 times the current, stel- 
lar tidal radii of the four dSphs. To gauge how important this pa- 
rameter is, we note that for Md = 10^, 10^° and 10"Afo, with 



V — 300 km we can double bmax and only induce a 10%, 20% 
and 45% respective increase in the Coulomb logarithm. Therefore, 
if one wishes to consider a larger bmax than we do here, simply sur- 
mise that the true Md is slightly lower than the used one. However, 
for realistic dSph masses « IO^'Mq) this effect is small. 



2.2 Milky Way mass model 

The dSphs mainly orbit in the halo of the Milky Way at Galac- 
tocentric radii >~ 20 kpc, therefore, the inner rotation curve - 
which is a hotly debated topic (see e.g. IBovv et al.ll2009l and ref- 
erences therein) - is of little concern. The main constraint we are 
interested in is the rotation curve measured at large radii from the 
Milky Way. The two article s we pay particula r attention to are those 
of G nedin etal.l <2010) and lXue et alj (I2OO8I). 

We use the same potentials as IXue et '^( 120081) for the bulge, 
disk and dark matter halo and we do not repeat them here for 
brevity. We choose Mdisk = 5 x 10^° M©, Mtuige = 1.5 x 
10"' M0, but like lMcGau"ghl ( I2OO8I) we choose the scale-length 
of the disk to be smaller (in this case b — 2.8 kpc, not 4 kpc 
as IXue et al]|2008l choose) in order to increase the rotation speed 
near th e solar radius to the expected Vc{rQ) « 220 kms~^ (see 
iMcGa ugh 2008). Of course we do not need it to be particularly ac- 
curate since we are using a spherical potential to model the disk 
and the fact is it does not influence the orbits of dSphs at 20 kpc, 
never mind beyond 100 kpc. 

Given that we are testing the ACDM model , it makes no 
sense to use a logarithmic halo, so we only use a iNavarro et al.l 
I1997I halo (hereafter NEW). We note, however, that a logarithmic 
halo has by definition stronger gravity in the outer parts of the halo 
and so the infalling satellites accelerate to larger speeds than with 
an NFW halo, which is counterproductive to the need to lose orbital 
energy, as we shall see later. 

To give the extremes of the permissible rotation curves, we 
take Fig 3 o f lGnedinetal] l l201(]|) which plots the allowed range of 
rotation curves between 20 and 80 kpc. From this we chose four 
halo models: one to represent a rotation curve that is roughly in the 
centre of the allowed range of velocities from 20-80 kpc (model 
I); one to represent the largest viable rotation curve (model II), an- 
other to represent the lowest viable rotation curve (model III) and 
a final one that falls from a large rotation speed at 20 kpc to a 
low rotation speed at 80 kpc (model IV). The NFW concentration 
parameters, c, and densities, ps, for the four models are given in 
Table 1 and we note that the virial radius at redshift zero, r2oo, is 
fixed to 275 kpc. Here r2oo refers to the radius at which the en- 
closed mass is two hundred times the critical density of the Uni- 



verse, A/200 ~ 200pc X iTrrloQ. We vary only the concentration 
parameter, c, for models I-III, but also the density, ps for model IV, 



('•200/c) 



. The rotation curves of the four models 



where p — ps t 7—r v2 

can be seen in Fig[T]as well as their dark matter densities and their 
typical DF decelerations for some reasonable parameters. 

Finally, increasing-mass Milky Way models are introduced in 
i]3.4l to account for the infall of dark matter onto the early Milky 
Way. 



2.3 dSph Proper Motions 

The crux of the whole issue revolves around the knowledge 
of the proper motions of the four dSphs, and of course their helio- 
cen tric radial vel ocities and distances from the Galactic centre (see 
e.g. lMatedlT998l) . The Galactocentric distances to the four dSphs 
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are given in Table |2] and all have relatively small errors, so we do 
not admit any variation in these current positions. 

The main uncertainty in accounting for the past orbits of the 
dSphs co mes from the pr oper motio ns. The proper mot ions of For- 



nax (Pia teket alJ|2007al) . Sculptor i Piatek et al 



aPiatek et ahfeoOSh and Carina Piatek et al 



2006 ). Ursa Mi- 
2003h . measured 



in those papers, were converted into Galactic rest frame tangential 
and radial velocities (by the original authors) and are given in Ta- 
ble[2] Clearly there is a significant amount of uncertainty associated 
with these velocities, and they have the potential to drastically alter 
the orbits. 



2.4 Orbit Evolution 

Our strategy, given the information available, is to elucidate 
the likelihood that the four dSphs fell into the Milky Way along a 
filament and ended up with their current positions and velocities. 
To do this, we begin from their current positions and velocities and 
trace back their orbits (accounting for DF) to see if, for a given 
dSph mass and Galactic potential, it is possible for the dSph to 
reach a significant distance from, or even escape, the Milky Way. 

The equations of motion are given in the appendix. 



3 RESULTS 

Our condition for a viable dSph candidate arriving by fila- 
mentary infall is the following: when we trace its orbit back in 
time, starting from its current position and range of possible veloc- 
ities, it must be able to reach a significant distance from the Milky 
Way within 12 Gyr. The maximum trace back time is curtailed to 
12 Gyr because any dSphs falling in before 12 Gyrs would have 
their positions with respect to the other dSphs randomised such that 
they could not account for the disk of satellites. This distance must 
be substantially greater than the current orbital apocentre found as- 
suming a fixed orbit with no changing Galactic potential and no 
dynamical friction. If this were satisfied, then it would suggest that 
dynamical friction is strong enough to decelerate an infalling the 
dSph and allow it to be consistent with its current proper motions. 

The variables discussed above include the Galactic potential, 
the proper motions and dSph mass. We have four Galactic poten- 
tials (parameters in TablefTJ, vary the proper motions between —3a 
and +3a of their mean measurements and consider dSph masses 
up to 10^^ M0. The aim is to find the minimum mass permissible, 
from constraints on the Galactic potential and proper motions, and 
then to judge the credibility of such a mass given the information 
we have on dSph masses from recent, independent observations. 

In Fig[2]we show the decay of a set of orbits for the Ursa Mi- 
nor dSph for various proper motions and using Galactic model I. It 
shows how critical the pericentre is for allowing the dSph to lose 
enough orbital energy to satisfy the current position and velocity 
after infall. Although the solid and dotted lines (mean measure- 
ment and mean-l-lcr error in proper motions) do not efficiently lose 
their orbital energy, the dashed line (mean subtracting the la error) 
does. Therefore, if Ursa Minor was observed to weigh upwards of 
5 X 10^" Mq and if the proper motions were confirmed to be la 
lower than the current mean measurement, then it would be a good 
candidate for filamentary infall. 



3.1 Dependence on the Galactic Model 

In Fig [3] we plot the largest apocentre reached within 12 Gyr 
against the dSph mass for the Ursa Minor dSph making use of all 
four Galactic potential models. The different linetypes, for Fig[3]- 
[S] correspond to the proper motions used: the solid line is the mean 
proper motion as defined in Table |2l the dotted and dashed lines 
correspond to izla and +3a/ — Xa errors respectively. There is a 
slight complication to this because in the cases of Carina and Ursa 
Minor, it makes no sense to consider proper motions as low as 3a 
below the mean measurement. The reason for this is that using the 
mean measurement minus la for Carina and Ursa Minor respec- 
tively, the pericentre is ~ 10 kpc. We stop at 10 kpc because as 
one can see from Fig|2l the pericentre does not vary considerably 
during a Hubble time, rather it is the apocentre that decays. There- 
fore, if a dSph had orbited to 10 kpc, we would expect it to be 
tidally destroyed. 

The closer a dSph approaches, the more efficient DF is, thus 
the easier it is to be captured. Interestingly, if the proper motions 
are found to be at the higher end of the error range, it is impossible 
- regardless of mass - to fall in and finish with the correct position 
and velocities. 

It is obvious from Fig[3]that there is a clear trend in all Galac- 
tic potentials for the largest apocentre to increase with dSph mass. 
This trend simply derives from the proportionality be tween the DF 
decele ration and the satellite mass (Eq[TJ - see also iTormen et al.l 
( Il998l) . Regardless of the Galactic potential, Ursa Minor would 
have required a mass greater than 10^'^ Mq, even at la below the 
mean, to have fallen in along a filament from more than 200 kpc. 



3.2 Dependence on Proper Motions 

Given that the Galactic model only weakly influences the abil- 
ity to lose orbital energy, we plot in Fig|4]the largest apocentre ver- 
sus dSph mass for each of the dSphs, but only for Galactic models 
I and II. 

Consider Fornax and Sculptor with Galactic model I (as well 
as Ursa Minor with Models I and III, FigO. For the +3(7 proper 
motions (where the proper motion error is added in parallel to the 
mean) the largest apocentre is constant with dSph mass. In this 
case, the dSph can always orbit beyond 700 kpc. We obtain a sim- 
ilar result for the +lcr where the largest apocentre is now between 
200 and 300 kpc. The largest apocentre also is almost indepen- 
dent of the dSph mass for Models III and IV, which are not shown, 
because these models have significantly lower circular velocities 
beyond 30 kpc (cf. Fig[T]( hence lower escape speeds. 

In contrast, it is the tendency of the majority of the cases 
where the proper motion is subtracted from the mean measure- 
ment in Fig|4]to increase their maximum apocentre with increasing 
dSph mass. In fact, the orbits referring to these lines begin 12 Gyr 
ago with large apocentres, but DF, whose efficiency increases with 
mass, reduces their orbital energy and they finish with much short- 
ened apocentres. This is a far more plausible scenario than the or- 
bits corresponding to the +2i-a lines. In this latter case, the dSph 
would still today be orbiting to beyond 700 kpc and this is likely 
to be inconsistent with the positions of the eight classical dSphs, 
which are observed within ~ 250 kpc of the Milky Way. In fact, 
out of the four we investigate here, three are within 90 kpc, with 
Fornax still only at 138 ± 8 kpc. Therefore, the probability for both 
Fornax and Sculptor to be on a long orbit to 700 kpc and being 
observed at 138 kpc and 87 kpc is small. Thus, it is reasonable 
to consider only scenarios where the dSphs today have relatively 
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small apocentres (say less than 200 kpc), hence excluding Galactic 
models III and IV. In addition, we can exclude Model II for Ursa 
Minor, Sculptor and Fornax (top-right panel of Fig[3]and right hand 
panels of Fig |4]l, because they have never been beyond 500 kpc 
from the Milky Way over the last 12 Gyr. 

It would appear, therefore, that there is no scenario which 
can reconcile Fornax with the filamentary infall, regardless of the 
proper motion, assuming it is within 3ct of the mean measurement. 

In the case of Sculptor, to have fallen in from far outside the 
virial radius, using the most favourable proper motion at Scr below 
the mean, the dSph would need to weigh more than IO^^Mq for 
Models I and II and at Icr it is barely possible regardless of dSph 
mass. 

Ursa Minor and Carina would both require masses of the or- 
der of 10"' A/0 and these masses come with the condition that the 
proper motions are at least la below the mean measurement. The 
reason for this, in contrast to Fornax, is that lowering the proper 
motions results in very low pericentres allowing much orbital en- 
ergy to be shed. 

A subtle point worth noting is why the largest apocentre is not 
a monotonically increasing function of dSph mass. This is linked 
to the number of pericentric orbits. By increasing the dSph mass 
slightly, the number of orbits can decrease by one, hence even 
though the DF strength at pericentre increases, the number of times 
it orbits to pericentre within 12 Gyr has dropped. 



3.3 Significance of dSpii mass 

In the previous section we confirmed that, using this simple 
model, there is very little chance that Fornax fell into the Milky 
Way's potential along a filament. Furthermore, as long as the proper 
motions of Sculptor, Carina and Ursa Minor are greater than or 
equal to the mean, it is very doubtful that their orbits could result 
from filamentary infall. However, if the proper motions co-conspire 
to be low at more than Icr, there is the possibility that dSphs with 
masses of the order of 10^" Mq could have fallen into the Milky 
Way's potential and ended up with their current positions and ve- 
locities. 

The problem with this scenario is that it is unlikely that the 
dark m atter halos of Sculptor, Ursa Minor and Carina are so mas- 
sive. In Angus ( 2010), cored dark matter hal os wer e fitted to line- 
of-sight velocity dispersions of Walker et'al] j2007l) . At 2 kpc the 
enclosed masses of Carina, Fornax and Sculptor were shown to be 
3, 3 and 2 x 10^ M q, Ursa Minor was not studied, but has a simi- 
lar mass profile (cf. lWalker etal]|2007h . The cored density profiles 
have the same asymptotic density profile as NFWs (i.e. oc r^^) and 
hence the cumulative ma sses at the virial radi i are not too different. 
In fact, the NFW fits from i Walker et alj ( l2007h have masses at 2 kpc 
that are approximately 0.8, 3.2 and 2x lO^M© respectively. 

Taking Carina, if we assume that it has a totally unphysical 
density profile such that the density stays constant at the 2 kpc 
value and integrate the total mass, we can find the minimum extent 
of the dark matter halo. Since the minimum dSph mass for filamen- 
tary infall is lO^'^Mo andp(2fcpc) = 2.8x lO'^A/ofepc"^, we find 
that Md = -^pr^ meaning r ~ 10 kpc. This is of course the best 
case scenario and it is more likely that the mass would need to be 
integrated beyond 20 kpc. Such diffuse halos would never survive 
the first pericentre passage and moreover, the DF formula (EqlTJ 
would be overestimating the effect for small impact parameters. 



3.4 Host Halo Accretion 

Given that we integrate backwards in time to 12 Gyrs ago, it 
is not satisfactory to exclude the accretion of mass to the host halo. 

To be able to elucidate the effect of this on our results, we 
followed the model of Wechsler et al. ( 2002) where the total mass 
of a cold dark matter halo is seen to grow with decreasing redshift. 
They suggest a prescription where the mass at any given redshift, 
M{z), is related to the mass M200 at redshift zero by M{z) — 
M2ooe~°'^ to which we impose a floor on the minimum mass, such 
that M{z) cannot drop below 5% of A/200 ■ Here a is the parameter 
that controls the rate of accretion and A/200 is the virial mass at 
redshift zero. 

If we recall that pc — ^^^q then we can equate the virial 

radius with the accretion rate by riooi^) ~ '^iooh'm^ ' 
the accretion, the gravity imposed by the Milky Way varies in pro- 
portion to both e~°^ and a factor coming from the growing virial 
radius with decreasing redshift. The parameters that control the DF 
- and the dark matter density p (see EqlTJ - vary as a conse- 
quence of this prescription and we also impose the same accretion 
rate on the stellar mass of the Milky Way (described in ^2.2\ . 

It is suggested that typical values of a for Milky Way sized 
halos are 0.5-0.7. We tested a fully in the range to 0.87 and plot 
in Fig[5]the results for a — 0.63 for Carina, Fornax, Sculptor and 
Ursa Minor, whilst limiting ourselves to Galactic model I. The left, 
middle and right hand panels show the maximum apocentre within 
5, 7 and 12 Gyrs respectively. 

At 5 Gyrs (see left hand panels of Fig|5j and within l-cr of 
the mean measurement of the proper motions, the results described 
in the previous sections remain valid, because there is only a slight 
difference between the maximum apocentres found here and those 
in Fig |4] that were derived considering a 12 Gyrs time interval and 
no accretion. Most apocentres are larger than those shown in Fig|4] 
because of the weakening potential with increasing redshift, rather 
than to the fact that the dSphs are escaping. Excluding Fornax, for 
masses smaller than 10^" Mq the dSphs are only orbiting up to 
250 kpc. 

The maximum apocentre of Fornax itself, at +l-a from the 
mean measured proper motion, does increase significantly by 
5 Gyrs. In Fig[5] left hand panel, one can see the +3-a line of For- 
nax reaches 400 kpc at low masses, whereas in Fig|4]the same line 
reaches only 250 kpc. By 7 Gyrs (middle panel of FigO this line 
has reached 500 kpc suggesting it could be in the throes of escap- 
ing. In addition, at 3-cr from the mean measurement, all four dSphs 
could already have escaped within 5 Gyrs, although this scenario 
would suffer the same flaws as stated in i i3.2l and H3.3\ either the 
dSphs are unlikely to be observed within ~ 140 kpc of the Milky 
Way when their orbits have apocentres larger than ^ 700 kpc or the 
required large mass of their dark matter halo is inconsistent with 
current mass estimates. sec:longorbs By 7 Gyrs (the middle panels 
of FigO the maximum apocentres have increased somewhat. Ursa 
Minor's maximum apocentres only increases very mildly, however 
the -1-1 -(7 line for Carina shoots up from 200kpc to 400kpc. This is 
a criteria for escape since its apocentre was stable at 200 kpc for 
5 Gyrs and between 5 and 7 Gyrs its apocentre increased by a fur- 
ther 200 kpc. Other dSphs that show apocentre increases of at least 
200 kpc are Sculptor (at -3-a and at the mean proper motion) and 
Carina (at +l-a from the mean proper motion), but certainly not 
Ursa Minor. The crucial point to notice here is that this only begins 
to happen after 5 Gyrs. 

Within 12 Gyrs (right hand panels of Fig|5} it is clear that 
the dSphs can reach a significant distance from the Milky Way. 
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One can see that by 12 Gyrs for the Carina dSph, it is possible to 
reach 1 Mpc from the Milky Way, even with a reasonable mass 
(10* M0). Likewise, Fornax and Sculptor can reach 800 kpc, but 
Ursa Minor can only reach 500 kpc for a reasonable combination 
of proper motion and mass. 

Considering that 12 Gyrs ago the virial radius of the Milky 
Way is only a fraction of its current value, it can be concluded that 
all the dSphs can infall, but not more recently than 5 Gyrs ago. 
Whether reaching 500 kpc or 1 Mpc from the Milky Way in 12 Gyrs 
corresponds to a large enough separation to be compatible with fil- 
amentary infall is not ascertainable with these simulations. Alter- 
natively, they may all have proper motions larger than the mean 
measurement by more than 1 — a and it could merely be a coinci- 
dence that these four dSphs are found within 140 kpc of the Milky 
Way whilst maintaining the long orbits to the outer halo they would 
have as a result of their large proper motions. The feasibility of the 
former scenario, namely an early infall, is not studied here, but has 
been extensively discussed and dismissed in iKlimentowski et al.l 
1 I2OIOI) , because the organized distribution of the dSphs we observe 
today is unlikely to survive such a long period of time. 



4 CONCLUSION AND DISCUSSION 

Here we took the positions and velocities of the four classi- 
cal dwarf spheroidal galaxies of the Milky Way that have mea- 
sured proper motions, and sought to discover if these final con- 
ditions were consistent wi t h the filamentary infall model pro - 
posed by 'Kang et al.| ( [200^ : iLibeskind et alj l l2005l l2007l |2009|) : 
IZentner et al. (2005J to offer a non-convoluted explanation of the 



sxpli 

so-ca l led disk of sate l lites (|K roupa et al] l2005l : iMetz & Kroupal 
I2OO7I ; iMetz et al.ll2009l : lLvnden-Bell,198j) . 

Using an analytic prescription of dynamical friction, and a 
simple leap-frogging time step updater with negative time steps, 
we were able to trace back the orbit for different masses of dwarf 
spheroidal galaxies in various Milky Way potentials. 

It was found that the proper motions and dwarf spheroidal 
masses are the two key parameters and that the Milky Way poten- 
tial model was rather subsidiary. The proper motions were varied 
up to ±3a from the mean measurement and all reasonable dwarf 
masses were considered. 

We argued that no conditions could be found such that For- 
nax might be consistent with the simple filamentary infall model. 
Moreover, the only scenario for Sculptor, Ursa Minor and Carina to 
be consistent with the infall hypothesis is for their proper motions 
to both be significantly lower than the mean, and for them to have 
total masses of more than 10^^ Mq. The problem with this scenario 
is that it is in stark conflict with the measured mass profiles of the 
dSph using the internal line-of-sight velocity dispersions of their 
stellar component. Realistically, the largest mass we could expect 
is 10^ Mq which is not conducive to the infall scenario. 

To further solidify our conclusion, we incorporated the effects 
of host halo accretion during the formation of the Milky Way's dark 
matter halo over the past 12 Gyrs. We found that this, by defini- 
tion, allowed the dSphs to reach larger apocentres and allowed all 
dSphs to move well beyond the Milky Way's virial radius (which 
decreases with increasing redshift). We also showed that all four 
dSphs remained bound at 5 Gyrs, unless the proper motions are all 
incorrect at more than \-o. Therefore, the filamentary infall sce- 
nario strongly depends on the dSphs falling in at redshifts larger 
than unity when the Milky Way is in its inf ancy. This scenario has 
been studied bv IKlimentowski et alj feOlOh amongst others and is 



not favoured due to the randomisation of the orbits of the many 
dSphs over the large intervening time. 

There is the possibility that a cold dark matter halo weigh- 
ing more than IO^'^Mq, that bound together several dwarf galaxies 
fell in, and in that way the dwarfs could lose enough orbital en- 
er gy. This however has t he same problem as the gr oups discussed 
in lMetz et ai] ( I2OO9I) and lKlimentowski et al.l j201(]|) . 

In addition to these points, we wish to point out certain logical 
disc repancies in some att empts to solve the Disk of Satellites prob- 
lem. iDeason et alj ( l201lh demonstrate that the anisotropic spatial 
distribution and angular momentum bias of the satellites at z=0 can 
be attributed to them being accreted relatively recently compared 
to the majority of the dark matter of the Milky W ay. Whereas, on 
the other hand, . Nichols & Bland-HawthomI ( I2OI 1ft claim to be able 
to account for the radial distribution of the dSphs if and only if they 
are accreted between z = 3 and z = 10. This example demon- 
strates that only small aspects of the whole problem can be solved. 
For example that the satellite galaxies of the Milky Way are gas 
poor can be understood if their gas was stripped, but this then re- 
quires them to have fallen in a long time ago. On the other hand, to 
explain the Disk of Satellites the dSphs must have fallen in recently. 

It would appear, therefore, that on the balance of this evidence, 
the filamentary infall model proposed to orchestrate the positions 
of the classical Milky Way satellites is untenable. We suggest that 
proper motions for the other four classical dwarf spheroidal galax- 
ies would greatly strengthen our conclusion, as would greater ac- 
curacy on the current proper motions. Any future proposition to ex- 
plain the disk of satellites in the concordance cosmological model 
must take care to adhere, not only to the highly organised positions 
of the satellites, but also to their proper motions. 

It is to be emphasised though that even in the presence of new 
proper motion measurements, it is not clear how the current conclu- 
sions can be changed. The new proper motion measurements would 
need to significantly deviate from the current best values, and the 
masses of the satellites would have to be made at least ten times 
greater than observed in order to give an infall solution. Both to- 
gether do not appear to be feasible, and so the conclusions reached 
here that the satellites cannot have fallen in appear to be robust. 

A further hurdle to the infall model is that while the gas may 
be distributed in a thin filament, the dark matter substructures, and 
notably those as heavy as 10^^ Mq have filament thicknesses com- 
parable to the diameters of the dark halo hosts they are feeding. 
Infall of dark matter sub-structures hosting satellites is therefore 
more likely to be nearly isotropic, or at best will lead to a tri-axial 
structure but with dimensions much thicker than the 10 — 30 kpc 
rms thickness of the disk of satellites. 

When this e vidence is added together with the elegant work 
of lPeebles & Nu sser (2010) that showed the strongly contradictory 
numbers of galaxies in the local void an d the con tinuing failure of 
cusps in low surface brightness galaxies jde Blokli2010.) , as well as 
the baryon i c TuU y-Fisher relation for gas rich galaxies described in 
iMcGaughl 1 I201TI) the deck is seriously stacked against the concor- 
dance cosmological model at the galactic scale. 

In certain m odified gravity theories such as M odified New- 



tonian Dynamics jMilgroml 



ISerra et alj|2010l : iKosowskv 



I995I a nd more recently lAngusll200^ : 



2OI0I) , where cold dark matter is re- 



placed by boosting gravity in regions of weak accelerations, the 
possibility of a simple origin of the dwarf spheroidals can hold be- 
cause they can (and must) be tidal dwarf ga l axies jX roupa 199^; 
Kroupa et alj I2OIOI: iBoum aud et al.l 12001 iGentile et aL ,20071: 
Milgrom 2 007 L lT^et & Com bes 2008) for which the correlated po- 
sitions and low orbital velocities are perfectly natural. 
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Figure 1. A series of important quantities are plotted against Galactocentric radius for tlie four different Milky Way mass models ( where the only differences 
are the dark matter halos). In panel (a) is the circular velocity out to 1 Mpc, (b) shows the dark matter density and panel (c) shows the deceleration due 
to dynamical friction for a dwarf satellite mass of 10^°Mq moving at 250 kms"\ The linestyles for models I, II, III, IV are solid, dotted, dashed and 
dot-dashed respectively. 





Figure 3. Here we plot the largest apocentre, as a function of mass, within 12 Gyrs, for the orbit of Ursa Minor in the four Galaxy models as defined in Table[T] 
The different linetypes correspond to the proper motions used: the solid line is the mean proper motion as defined in Table |2] the dotted and dashed lines 
correspond to ztltr and +3(t errors respectively. Regardless of the Galactic potential, Ursa Minor would have required a mass greater than IO^^Mq, even at 
1(T to have fallen in along a filament. 
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Figure 4. Here we plot the largest apocentre, within 12 Gyrs, for the orbits of Sculptor, Fornax and Carina with Galaxy models I and II - as defined in Table[T] 
The different Hnetypes coirespond to the proper motions used: the solid line is the mean proper motion as defined in Table [2] the dotted and dashed lines 
conespond to ztlcr and ±3(T eiTors. Caiina has only one dashed line as explained in i|3.2l 



Model 




c 


I 


8 


13.5 


II 


8 


12.0 


III 


8 


16.5 


IV 


340 


75.0 



Table 1. In the left hand part of the table we present the scale densities and 
concentration parameters for the four Galactic dark halos. r200 is always 
set to 275 kpc at redshift zero. 



dSph 


ro [kpc\ 




Vyn (kms -*-) 


Fornax 


138±8 


-31.8±1.7 


196±29 


Sculptor 


87±4 


79±6 


198±50 


Ursa Minor 


76±4 


-75±44 


144±50 


Carina 


101±5 


113±52 


46±54 



Table 2. Here we give the Galactocentric distances and velocities of the 
dwarf spheroidal galaxies. For Fornax, Sculptor and Ursa Minor, our Vxq 
con'esponds to Piatek et al's Vr an d our Vy„ to t heir V t. For Carina the 
proper motion comes directly from IPasetto et al] feOllh . Distances come 
from ,Matea I .199&) . 
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Figure 5. Like for Fig|4]except we include tlie effects of liost halo accretion (see i|3.4t . All panels use a = 0.63, but the left, middle and right hand panels are 
the largest apocentre after 5, 7 and 12 Gyrs respectively. We only plot Galactic model I, however, the results are the same for model II. 
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APPENDIX A: INTEGRATION SCHEME 

In order to be clear about the way the dynamical friction af- 
fects the equations of motion, we write them in full here. 

= Xi+Va:-At/2 

y^+l/2 = yi + Vy^At/2 

Xi+l = a;i + i/2 + ^a;, + i At/2 
Vi+l = yi+l/2 + Vyi+iAt/2 

(A2) 

where rt = -y/a;? + yf, Vi = ^/V^i~-^-V^, and At is the time step 
which is negative. 

Note that initially we assume all dSphs are at s = ro and y = 
0, where ro is the current Galactocentric distance and therefore that 

is the current (and initial) radial component of velocity and Vyg 
is the tangential component, as given in Table |2] The calculations 
are done in the plane of the orbit just to minimise the degrees of 
freedom, and does not affect the results in any way. 



